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finite isospin density. 

PACS numbers: 12.39.Fe; 11.30.RD; 12.38.Lg; 



E-mail: zhaozhang@pku.edu.cn 
2 E-mail: yxliu@pku.edu.cn 



1 



1. Introduction 

The phase structure of QCD at non-zero temperature and baryon chemical potential has been 
intensively investigated throughout the last decade. In reality, dense baryonic matter obeys 
an isospin asymmetry i.e., in the case of two light flavors, the densities of u and d quark are 
different. In order that QCD adequately describe the isotopically asymmetric matter, such as 
compact star, isospin asymmetric nucleon matter and heavy ion collisions , usually the isospin 
chemical potential fii = (jj, u — A*d) is introduced in the theory[TJ[2]. Different approaches, such as 
Lattice QCD[3ll3], chiral perturbation theory [TJ [21 Ladder QCDp], Nambu-Jona-Lasinio 
type model [H [91 [THl [HI [E], and random matrix model[13] has been used to explore the QCD 
phase structure at finite isospin density. It has been widely confirmed that there is a phase 
transition from the normal phase to the pion superfluidity phase which is characterized by a 
pion condensate (d'j^u + H.c.) at high enough isospin chemical potential. It is also found that 
the kaon superfluidity phase characterized by kaon condensate (sj^u + H.c.) appears at high 
isospin and strangeness chemical potential in the three light flavors case [21 [TT], 

The previous studies on the effects of finite isospin chemical potential and strangeness are 
mostly focused on two types of condensates {uj^d + H.c.) and (575 d + H.c), which are order 
parameters for the corresponding superfluidity phase transitions . It is generally believed that 
the vacuum of QCD has complicated structure and it is expected that all gauge invariant 
Lorentz singlet local operators built of the quarks and/or gluons have non- vanishing vacuum 
expectation values according to QCD sum rules |141 I15j . For example, the well-known low- 
dimensional condensates, such as quark condensate (qq), gluon condensate g 2 (GG) , mixed 
quark gluon condensate g(qcrGq) and four quark condensate (q^iqq^2q), play significant roles 
in the hadronic studies based on QCD sum rules. 

Due to the presence of the flavor mixed condensates (uj^d + H.c.) and {s^d + H.c.) at finite 
isospin and strangeness chemical potential, it is natural to expect that there should exist other 
new types of flavor mixed condensates induced by isospin chemical potential and strangeness 
chemical potential according to Operator Product Expansion (OPE). Besides the pion conden- 
sate and kaon condensate, the possible low-dimensional flavor mixed condensates are mixed 
quark-gluon condensate g(dj^aGu + H.c.) induced by isospin density, and g^sj^aGu + H.c.) 
induced by strangeness density (for convenience, we call the former pion mixed quark-gluon 
condensate and the later kaon mixed quark-gluon condensate). In addition, new forms of four 
quark condensates, such as (qj5TiTiqqT2q), may also appear in OPE. We expect these induced 
low-dimensional condensates also play important roles on the hadronic physical observables in 
the framework of QCD sum rules. 

It is well-known that, in the chiral limit, both chiral condensate and mixed quark-gluon 
condensate are ideal order parameters for the chiral phase transition of QCD. Similarly, pion 
mixed quark-gluon condensate and kaon mixed quark-gluon condensate can play the roles of 
order parameters for the pion superfluidity phase transition and the kaon superfluidity phase 
transition at finite isospin chemical potential and strangeness chemical potential, respectively. 
Though both pion condensate and pion mixed quark-gluon condensate can be used as order 
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parameters to describe the pion superfluidity phase transition, they reflect different aspects of 
the non-perturbative structure of the ground state: the former reflects the correlation between 
different flavors with color-singlet component, while the later reflects the correlation between 
different flavors with color-octet components. Therefore, pion mixed quark-gluon condensate 
will give new and important information on the pion superfluidity phase transition. The same 
thing is true for kaon condensate and kaon mixed quark-gluon condensate. 

Therefore, it is interesting to investigate the thermal and dense properties of these new 
types of low-dimensional condensates and their effects on the physical hadronic observables. 
Since the Global Color Model (GCM)PH El El [TO] is an effective Q uark and g luon fields 
theory and has been successfully used to investigate the property of the traditional mixed 
quark-gluon condensate [20l [21] and other QCD condensates [22], we will adopt this model to 
explore the thermal and dense properties of above induced mixed condensates. In this paper, 
we only consider the possible pion superfluidity phase transition and pion mixed quark-gluon 
condensate. 

2. Mean field theory of GCM at finite isospin chemical potential 



In Euclidean metric, with {7^,7^} = 25^ and 7+ = 7^, the generating functional of GCM 
with quark and gluon degrees is 

Z[J,rj, V ]= J DqDqDAexp(-S GC M[q, ?, A% +rjq + QV + (1) 



with the action 



S G CM[q,q,A<$ = J (qfr.d + M-igA«^)q + ±A«D-i(id)AZ). (2) 

The essence of GCM is that it models the QCD local gluonic action J F^ V F^ V which has local 
color symmetry, by a highly nonlocal action which has a global color symmetry. The main 
aspects of GCM have been reviewed in [T7 1 fT8 ] fl9]. 

Integrating over the gluon degrees, the partition function of the GCM with two quark flavors 
at finite baryon and isospin chemical potential(we only consider the case with temperature 
T = in this paper) has the form 

Z(ji, m) = j Dq(x)Dq(x)exp[- J q{x)[y^d^ + M - ^74 - SfiT 3 j4]q(x) 

-\ JJj*(x)g 2 D^(x - y; (3) 

where M = diag(m u ,md), J x = J d x and j®(x) = q(x) r y fJ- -^-q(x). In Eq.®, Ti(i = 1,2,3) 
are the Pauli matrixes in flavor place, [i = /xb/3 is the chemical potential associated with 
baryon number and the quantity b[i is a half of the isospin chemical potential, i.e. b[i = m/2 
(In this paper, we only consider /ij > 0). The effective gluon propagator g 2 D llu (x — y; fi, fij) 
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is generally a (/x, /independent function, which is parameterized to model the low energy 
dynamics of QCD. 

In this study, we will take m u = m d = m. Evidently, the above Lagrangian is invariant 
under the baryon £/b(1) symmetry and the parity symmetry transformation P. In the case 
with fij ^ 0, the traditional isospin SUi(2) symmetry is reduced to Ui 3 (l) symmetry. Usually, 
the quark condensate (qq) is responsible for the chiral symmetry breaking of the ground state 
and does not spoil the parity and isotopical symmetry, while the nonzero pion condensate 
(qj5Tiq) breaks both the parity and isotopical symmetry of the ground state. 

Within the GCM formalism, the ground state of QCD is defined by the saddle point of 
the action and the quark gap equation at the mean field level is determined by the rainbow 
truncated quark Dyson- Schwinger equation (DSE) (The application of the DSE model to finite 
temperature and density is reviewed in |23j) 

S(p) = \J J^9 2 D,„(P ~ q)l»S{q) lu . (4) 
At finite (h, Hi) with u, d quarks, the inverse of quark propagator can be written in the form 



o-i/ n c -i/ \ . ( E uu (p, m) XwO, fi, Hi) 

\ Ti du {jp, h, m) ^dd{P, H, Hi) 



(5) 



where 

e-i/ \ f il -p + ilAWu + m \ , . 

S (p,/i,M/) = . > ( 6 ) 

\ ^7 • p + ^74^fc^ + m J 

with w u = (p4 + ifi + iSfj), Wd = (P4 + ifi — iSfi), 

Z aa (p, fj,, hi) = H ■ pA a (p,w u ,w d ) +ij4W a B a (p,w u ,w d ) + C a (p,w u ,Wd), (7) 
^ud{p,H,Vi) = ^ du (p, h, Hi) = i-y 5 D(p,w u ,w d ), (8) 

and A a ,B a ,C a ,D are momentum-dependent scalar functions. Nonzero C a and D are responsi- 
ble for the dynamical chiral symmetry breaking and isotopical symmetry breaking, respectively. 

Note that the possible diquark condensation is not considered here and only single Lorentz 
structure is concerned in H ud and T, du . The four matrix elements of the momentum dependent 
quark propagator 

S(p fj, hi) = ( ^(P'^Vi) S ud (p,H,Hi) ) / g x 
\ SduiP, fJ>, Ml) S dd {p,H,Hi) J 

take the form 

Sun = [(X d C u + D 2 C d ) - il ■ p(X d A u + D 2 A d ) - i lA (X d w u B u + D 2 w d B d )]/H, 
Sdd = [(XuCd + D 2 C U ) - i7 • p(X u A d + D 2 A U ) - m(X u w d B d + D 2 w u B u )]/H, 



Sud = ~iD[ l5 Y - , 75 7 • pS™ - i 7 574^r + 757 • Pl^^/H, ^ 
Sdu = ~iD[ 75 Y - z 7 57 • PSIT ~ ^Sjf 4 - 757 ■ FT^f "»]/H 
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with 





17574 _. 

ud 



(11) 



<75774 _ 
ad 



and 



X a = A 2 a f + B 2 a w 2 a + C 2 , Y = C u C d + A u A d f + w u w d B u B d + D 2 , 
H = X u X d + D 4 + 2D 2 (C u C d + A u A d f + w u w d B u B d ). 



(12) 
(13) 



With above decomposition, the gap equation can be expressed as 
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5 D^uip- qhnSij(q)j v . 



(14) 



Since each term of Eq. fllip is nonzero at finite /i/, more Lorentz structures with new scalar 
functions should also be considered in Y, ud and T, du to guarantee the self-consistent treatment of 
the gap equation. However, introducing more Lorentz structures will complicate the resolving 
of the gap equations. Just as the Lorentz tensor structure is not concerned in S aa at the 
traditional treatment of DSE at finite (T, //)[2$j, we suppose that i^D is the leading order 
term of Ti ud and T, du and other Lorentz structures has small impact on the determination of 
quark self energy. Because there are no structures 757, 7574 and 75774 in T, ud and T* du , the 
corresponding structures associated with Eq. ([TT]) in S ud and S du are ignored in the following. 
At least, this is a good approximation in the case with small fij. 

Due to the phenomenological nature of this effective theory, for simplicity, the Feynman-like 
gauge g 2 D^ u (p—q) = 5^ u g 2 D(p—q) was adopted in our calculation. With above approximation, 
the gap equation (fT4"j) is reduced to seven coupled integral equations, which is still complicated 
to solve. To get a qualitative understanding of the phase diagram and the structure of the 
ground state at finite a pedagogical model first introduced by Munczek and Nemirovsky [23] 
for the modelling of confinement in QCD is favored in this study. Munczek-Nemirovsky(MN) 
model has been extensively used to explore the properties of strong QCD both at zero (T, /i) 
and nonzero (T, ,u)f25ll2T]. which can always give qualitatively consistent results with the more 
sophisticated models. The effective gluon propagator of MN model takes the form 



with the single parameter 77 determined by tt and p masses in vacuum. The scale parameter 
r\ has relation with the string tension of QCD, and in the more real world, it should been 
the function of T and p. Using Eq.([T5|), the complete expressions of Eq. (fT4"|) are simplified as 



9 2 D, v {p -q)= 5*(2k)S 2 5\p ~ q) 



(15) 
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seven-coupled algebraic equations 

{A u - 1) = \r) 2 [X d A u + D 2 A d ]/H, (A d - 1) = \r] 2 \X u A d + D 2 A U ]/H, (16) 

(B u - 1) = \tf\X d B u + —D 2 B d ]/H, (B d - 1) = \if[X u B d + —D 2 B U ]/H, (17) 

C u -m = r, 2 [X d C u + Z) 2 Cy/tf, C d - m = rj 2 [X u C d + D 2 C U }/H, (18) 
= r) 2 D[C u C d + + w u w d B u B d + D 2 ]/F. (19) 

Eq. (|19p illustrates that there are two distinctive solutions to D: one characterized by D = 0, 
which describes the normal phase; the alternative, characterized by D ^ 0, which describes the 
pion superfluidity phase. The phase with small free energy is favored in nature. 

3. Thermal Potential and Condensates 

In GCM/DSE formalism, whether the normal phase or the pion superfluidity phase is stable 
is determined by evaluating the (/i, /i/) -dependent pressure difference 

SPfam) = P[v,fii,S[D + 0]] - P[>,mj,S[Z> = 0]], (20) 
where the pressure is calculated by using a steepest-descent approximation [26J : 



/' •/•./'•/'/.*: = -n[S] = //•/./, is ; - / , /.„ >:.s'. (21) 

Using the technique 

Det f C Z> ) = Det ( A ) Det ( B ) Det ( C ) Det ( C ' lDB ' 1 ~ ( 22 ) 
the pressure at finite (zx, zxj) can be expressed as 



P[S] = [ 2Ln[H] + 2N C f \D 2 [f(A u + A d ) + w u w d {B u + B d )\ + X u [p 2 A d + ^BJ 
+ X d [p 2 >l u + x/; 2 ^] + m[X d C u + X u C d + D 2 (C U + C d )]l 



(23) 



where f p = J ^J)i • Note that a constant term has been ignored in Eq. (|23p . Though the 
pressure (or thermal potential) calculated through Eq. (|23p is ultraviolet divergent, the pressure 
difference or the "bag constant" 8P |16j is finite. 

From the GCM generating functional, it is straightforward to calculate the vacuum expec- 
tation value(VEV) of any quark operator with the forms 

°n - Gi.AgUi^AgU) • • • GkA&ftJ, (24) 

in the mean field vacuum. Here AW stands for an operator in Dirac, flavor, and color space. 
The VEV of the operator O n has the form [27] 

{O n ) = (-ir£(-HA« ■■■^±S iljp(1) ...S injpin) ], (25) 



6 



where p stands for a permutation of the n indices. Based on formula (|25|) . the low-dimensional 
condensates, such as chiral condensate and pion condensate can been expressed as 



(uu) = -Tr DtC [S uu (x,x) - a uv (x,x)} = -N c Tr D [S uu (p) - a uv (p)], (26) 

J P 

(dd) = -Tr Di c[S d d(x,x) - a uv (x,x)] = -N c Tr D [S dd (p) - a uv (p)], (27) 

Jp 

(qij 5 Tiq) = -Tr D ,c,F[ij5nS(x,x)} = -N c Tr D [ij 5 S ud (p) + i^S du {p)]- (28) 



To get a convergent condensate integral, a subtracting term (Juv{p) which simulate the ultra- 
violet behavior of the quark propagator is introduced in the definition of the quark condensate. 
In the case with nonzero current quark mass and zero chemical potential, for MN model one 
has, with s = p 2 , 

A(s) = B(s) = l + ^, C(s) = m(l + -), (29) 
2s s 

with corrections of high order in (1/s). From these approximate expressions one can construct 

auv{3) = sA(8)* + CW (30) 

Since cruv(s) ~~ > m / s as s — > cxd and auv( s ) as s -> 0, this prescription will provide 
an absolutely convergent result with no need for a cutoff. The above definition can be gen- 
eralized to the case with nonzero quark chemical potential. There is no need to introduce a 
subtracting term in the definition of the pion condensate since D(p) keeps zero in the large 
momentum region in MN model. The formula for evaluating the 6-dimensional four-quark 
condensates can also be directly derived from Eq. (l25p . which is consistent with the vacuum 
saturation approximation at zero quark chemical potential. At finite isospin density, the new 
type four-quark condensates, such as (qi'JsTiqqq) , will appear in OPE. It is easily to prove that 
{qi^Tiqqq) ~ (qi'ysTiq} {qq) with the approximation that only 75 structure is hold in ^ u d(du) 
and S ud{du) , with (qq) = (uu + dd). 

Since the functional integration over the gluon field A^ is quadratic in the framework of GCM, 
one can perform the integration over gluon field analytically. Using the technique introduced 
by MeissnerpO], through the following integral formulae 

J VAe^ AD ~ lA+jA = e^i, 
I VAAe-^ AD ~ XA +i A = (jD)e^ D i, 

J (31) 
J VAA 2 e-^ AD ' lA+ ^ A = {D+ (jD) 2 ] e y D i, 
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the gluon fields vacuum average can be replaced by the quark current j° with the effective 
gluon propagator D(x — y). At the mean field level, according to Eq. (l25p . one can in principle 
obtain the VEVs for any gluon fields. This technique provides a feasible way to calculate the 
VEVs of operators with low-dimensional gluon fields such as the traditional mixed quark-gluon 
condensate and the isospin density induced pion mixed quark-gluon condensate. Since the 
number of terms produced by Eq. (|25p will increased rapidly with the number of gluonic fields, 
this technique is not suitable for the evaluation of the VEV of the operator involving high 
powers of gluonic field A. For instance for the gluon condensate (GG), which contains a A 
term, the calculation gets already rather involved. 

Applying the method described above, we obtain the expression 

giqi^na^G^Yl) = ~ 2iN c J ^[ d v9 2 D(y - x)]Tr DjF [S(y - x)i^T 1 a jjLV S(x - y)^ v \ 

+ l/.V, / / g 2 D(y - x)g 2 D(z - x)Tr DtF [S(z - x^na^ ^ 



ly J z 

xS{x- y)^S(y - z)7„]. 

The similar expression for evaluating the traditional quark-gluon condensate g(qa a uG a au -fq) 
can be got by replacing the structure i^T\(T au in (|32l) with a^ v . 
Using the Eq, (|14p and the formulae 



\a \b \c \a \c \b 

^c\. _ A c A c A c 

2 2 2 



-J° 



abc fabc 



r c f 



N c 5 aa , 



2 2 2 

the expression for the pion mixed quark-gluon condensate can be simplified as 

K 

2 



giqiKna^G^^-q) = h + 1 2 + h + h + h + h + h, 



where 



72 / §y 



(A u + A d - 2)]f + (B u - l)w 2 u + (B d - l)w 2 d 



X d A u + D 2 A d + X u A d + D A u 



X d B u w 2 u + D 2 w u w d B d + X u B d w\ + D 2 w u w d B u 



81 



h = — =Y 



D, 



2 . L H 



D 1 



(C u - m)(C d - m) 



2 .L H 



(C u - m){X d C u + D 2 C d ) + (C d - m)(X u C d + D 2 C U 
(A u - l)(X d A u + D 2 A d ) + (A d - l){X u A d + D 2 A u )]f, 



2 .L H 



(B u - l)(X d B u w 2 u + D 2 B d w u w d ) + (B d - l)(X u B d w 2 d + D 2 B u w u w d ) 



(33) 

(34) 

(35) 
(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
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Note that the scalar functions A a , B a ,C a , D are all momentum dependent in the DSEs formal- 
ism. The above integral expressions explicitly show that the pion mixed quark-gluon condensate 
is a order parameter for the pion superfluidity phase transition. By the way, it should be men- 
tioned that the expression in Eq. (l34p for pion mixed quark-gluon condensate is only valid in 
Feynman-like gauge. 

4. Numerical Results and Discussions 

To get a qualitative understanding of the finite [ii effects on the ground state, the numerical 
study below are all based on the MN model. Due to the defect of the supposed Lorentz 
structure of quark self-energy, the isospin chemical potential concerned below is limited within 
the range |/i//2| < 0.2(GeV). In this paper, for simplicity, we only concern the case with zero 
temperature and baryon chemical potential. 

A. The critical isospin chemical potential for pion condensate 

The effective field theory arguments [1] indicate that critical isospin chemical potential fi'j for 
the pion superfluidity phase is exactly the vacuum pion mass m n at T = [j,b = 0. Whether the 
MN model can algebraically reproduce this result within above formalism is discussed below. 

For the solution with Dj^O , the gap equation (fTTJj) is reduced to 

H = t] 2 [C u C d + A u A d f + B u B d w u w d + D 2 } . (42) 

In the neighborhood of fij c , one can probably neglect the D 2 and D 4 in above expression. 
According to gap equations (fl"6j) and (fTT|) . A u ^ and B u ^ are identical for D = 0. Therefore, 
the gap equation (|4"2|) can be further reduced to 

pi pi pi pi pi pi 

C(P + —)C{p - —) + A(p + —)A(p - — )(p + — ) • (p - —) 



X(p + ^-)X(p-^-) = r ] 2 



(43) 

where P' = (0,i/j,j) and the flavor subscript has been ignored. Within the feynman-like gauge 
and only considering the 75 structure, the Bethe-Salpeter equation (BSE) for the vacuum 
pseudoscalar amplitude Ti(p,P) in GCM takes the form [T7] 

2 r d 4 a P P 

r wj (p,P) = --J ^D(p-q)Tr D ,cA^r j + S(q+-)i l5 r j S(q--)]T 7rj (q,P), (44) 

where P 2 = —m 2 . Using Trp[T^~Tj] = 2, the BSE (|44p is greatly simplified in the MN model 



X(p+^)X(p 



P 



) = v 2 



A(p + ^)A(p - ^){p + |) • (p - ^) + C(p + ^)C(p - £ 



(45) 



which has the same form as the gap equation (|43p . Therefore, if the BSE (|44p can produce 

-P' 2 = m 2 and get the conclusion ^ = m n at 



the vacuum pion mass, we algebraically prove 
T = fi = . 
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However, the above proof is only true for the chiral limit case since the MN model doesn't 
support the pion bound state beyond the chiral limit if only the 75 structure is considered 
in the pseudoscalar meson Bethe-Salpeter amplitude [21] . Therefore, to produce the result 
\i\ = m-x beyond chiral limit in the MN model, other Dirac amplitudes beyond 75 -D should also 
be included in the off-diagonal term of the inverse quark propagator, which will make solving 
the gap equation more involved. 

Note that for other improved effective gluon propagators such as those used in [201 1211 [22] , 
the pion Bethe-Salpeter amplitude with only 75 structure is a good approximation to obtain 
the vacuum pion mass. In principle, using these improved effective gluon propagator to explore 
the pion superfluidity within the DSE formalism should give more quantitatively reasonable 
results in contrast with the simple MN model. However, the set of seven coupled algebraic gap 
equations (fT6j) - (Tl9j) will be replaced by a set of seven coupled integral equations which leads to 
the numerical calculation more difficult. 

For simplicity and to get a qualitatively understanding of the pion superfluidity within the 
DSE formalism, we still use the MN model in the following and only contains the 75 structure 
in the off-diagonal term of the inverse quark propagator. 

B. The Chiral limit 

T=0, a*b=0 

2 
1.5 



1 



0.5 


0.05 0.1 0.15 0.2 

W/2(GeV) 

Figure 1: The pressure difference 5P, pion condensate ir and pion mixed quark-gluon condensate 
Mix in the chiral limit. Bo = (O.lry) 4 , (qq)o = (uu + dd)o and Mixo = (310MeV) 5 are the bag 
constant, chiral condensate and mixed quark-gluon condensate of the vacuum obtained from 
the MN model with m = and r/ = 1.06GeV, respectively. 

In the chiral limit, there are four possible solutions according to the gap equations (|16|) - (|19|) : 

C u = C d = 0,D = 0; C u ^0,C d ^0,D = 0; 
C u = C d = 0,D^0; C u ^0,C d ^0,D^0, 
which characterize four possible phases of QCD at finite isospin density, respectively. However, 
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the solution with both nonzero quark condensate and nonzero pion condensate is not found 
in our numerical study. By comparing the corresponding free energies of the former three 
possible phases, it is found that the pion superfluidity phase is favored in the chiral limit at 
nonzero isospin chemical potential. It seems that this is a universal result which has been 
confirmed by many former studies [1, 2j. In this case, chiral symmetry is not broken and 
both quark condensate and mixed quark-gluon condensate disappear in OPE. In contrast with 
the vanishing of above chiral condensates, new condensates induced by isospin density, such as 
pion condensate and pion mixed quark-gluon condensate appear. Note that the new four-quark 
condensates, such as (qi'JbTiqqq) , also vanish in this case because these condensates factorize 
at mean field level within our formalism with the approximation used in Sec. 2. 



T=0,/i B =0 




W/2(GeV) 

Figure 2: The isospin chemical potential dependence of the ratio between pion mixed quark- 
gluon condensate and pion condensate from MN model with m = and rj = 1.06GeV. 

The m dependence of the pressure difference between the pion superfluid phase and the 
normal phase, pion condensate and pion mixed quark-gluon condensate is shown in Fig.l. 
For the case with \xi = 0, the zero pressure difference of 5P suggests that the nonzero pion 
condensate and nonzero quark condensate correspond to equivalent but distinct vacua, which 
is guaranteed by the chiral symmetry(A small quark mass will destabilize the superfluidity 
phase). In contrast with this event, for the whole domain of nonzero fij concerned, the positive 
pressure difference SP suggests the pion superfluidity phase is the stable ground state in the 
chiral limit for two-flavor QCD. Fig.l shows that the magnitudes of both induced condensates 
and the pressure difference are monotonically increasing functions of ///. 

Fig. 2 shows the iij-dependent behavior of the ratio of pion mixed quark-gluon condensate to 
pion condensate from MN model. In Fig. 2, The ratio ranges from 1.56(GeV) 2 to 1.66(GeV) 2 , 
which suggests that the induced mixed quark-gluon condensate has the same magnitude of the 
traditional mixed condensate in the concerning domain of /ij (in MN model, the ratio of the 
mixed quark-gluon condensate and the chiral condensate is 1.92 in the vacuum [2T]). Note that 
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numerical study suggests that the pion mixed quark-gluon condensate defined in Eq. (l32p is 
positive, which is consistent with the sign of the traditional mixed quark-gluon condensate. 

C. Finite current mass 

T=0, Hb=0 




W/2(GeV) 

Figure 3: The pressure difference SP ( plotted only in the region > /j,f ), quark condensate 
(uu), pion condensate tt and pion mixed quark-gluon condensate Mix obtained from MN model 
with m = 12MeV and 77 = 1.06GeV. 

With finite current quark mass m = 12(MeV)|24j. apparently, only two types of solutions 
exist: the normal phase with nonzero quark condensate and zero pion condensate and the pion 
superfluidity phase with nonzero pion condensate and nonzero quark condensate. In contrast 
to the chiral limit, for finite current mass, there is no solution with zero quark condensate to 
the gap equations due to the explicit chiral symmetry broken term in the lagrangian. 

According to the numerical study, in the small isospin chemical potential region fij < 64MeV, 
only the normal phase solution exists. For the region jU/ > 64MeV, there also exist the solution 
corresponding the pion superfluidity phase beside the normal phase solution. The stable ground 
state in the region jit/ > 64MeV is again determined by the difference of the pressure 6P in 
Eq.([23D. 

It is shown in Fig. 3 that at the point = 64MeV, the scaled pressure difference of two 
solutions is close to zero and for the region > 64MeV, the pressure difference is positive 
and monotonically increase with the increase isospin chemical potential, which suggests that 
the pion superfluidity phase is favored in > 64MeV region. From the ///-dependent behav- 
ior of both pressure difference 5P and the two order parameters, pion condensate and pion 
mixed quark-gluon condensate, one can judge the phase transition from normal phase to pion 
superfluidity phase is second order. This conclusion is consistent with the result obtained from 
lattice simulation and other model studies. 

Note that at zero fi, the scalar functions A(B,C) U in S uu and A(B,C)d in Sm at the same 
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Figure 4: The ///-dependence of the ratio between pion mixed quark-gluon condensate and 
pion condensate obtained from MN model with m = 12MeV and r/ = 1.06GeV. 



point (p,w u ,Wd) are complex conjugates, therefore the relation (uu)=(dd) always holds. Fig. 3 
shows that the magnitude of quark condensate monotonically increases with the increase of /// 
in the normal phase region, which is consistent with the dependence of the chiral condensate on 
the baryon chemical potential obtained within DSE formalism [23]. In the pion superfluidity 
phase, the magnitude of quark condensate monotonically decreases with the increase of /jlj, 
which is an anticipated result due to the monotonically increasing behavior of pion condensate 
with respect to fij. This behavior is similar to the dependence of the chiral condensate on 
the baryon chemical potential for the appearance of diquark condensate |23j . It is expected 
that the traditional mixed quark-gluon condensates, g(uaGu) and g(daGd), have the simi- 
lar ///-dependent behavior. Fig. 3 manifests a competitive relationship between the induced 
condensates and their corresponding traditional partners. 

In Fig. 4, we display the numerical result of the isospin chemical potential dependence of the 
ratio of the pion mixed quark-gluon condensate to pion condensate. Such a value ranges from 
2.2(GeV) 2 to 1.7(GeV) 2 , which is also close to the ratio of the traditional mixed quark-gluon 
condensate to quark condensate obtained in the vacuum [2 1J. The large magnitude of the ratio 
suggests that the induced mixed quark-gluon condensate is an important parameter within the 
QCD sum rules at finite ///. In addition, in contrast with the chiral limit case, the nonzero 
pion condensate and nonzero quark condensate suggests that the new four-quark condensates 
also have nonzero value in the superfluidity phase, even at the mean field level. 

The critical chemical potential 64MeV is relatively small in contrast with the vacuum pion 
mass. As mentioned above, the main reason for this discrepancy arises from the fact that only 
75 structure is considered in the off-diagram part of the inverse quark propagator, while the 
MN model doesn't support the pseudoscalar bound state beyond the chiral limit when only 75 
structure is contained in the pseudoscalar meson Bethe-Salpeter amplitude. One can expects 
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Figure 5: The dependence of the critical isospin potential ^ on the current quark mass from 
MN model with rj = 1.06GeV. The current quark mass ranges from to 20MeV. 

that this discrepancy will become small when we either adopt the other improved effective 
gluon propagators in the calculation or include other allowed Dirac structures such as 757 • p 
in the off-diagonal part of the inverse quark propagator within the MN model. 

The dependence of the critical isospin potential fij on the current quark mass is plotted in 
Fig. 5. It is shown that \x\ monotonically increases with the value of current quark mass. 
Though nj obtained in MN model with only 75 structure considered is markedly smaller than 
m n , the critical point is still roughly proportional to the square root of current quark mass in 
the range (0-20) MeV. 

5. Summary and Remarks 

Using a pedagogical confining model within the framework of GCM, we have qualitatively 
investigated the phase structure and condensates of two-flavor QCD at finite isospin density 
with zero temperature and baryon chemical potential. By solving the quark gap equation 
through the DSE formalism, we obtained that the truncated DSE type model supports the 
pion superfluidity phase transition at high enough isospin chemical potential. In contrast with 
the previous model studies, the obtained gaps responsible for both the chiral condensate and 
pion condensate are all momentum dependent within the DSE formalism, which are more close 
to the real world. In addition, some new types of low-dimensional condensates of QCD induced 
by finite isospin chemical potential, such as pion mixed quark-gluon condensate and mixed four 
quark condensate, are proposed and investigated in this paper. 

In the chiral limit with finite isospin chemical potential, the normal phase is unfavored and 
the pion superfluidity phase is the stable ground state. For the case with finite current quark 
mass, only the solution corresponding to the normal phase is found in the gap equations in 
the region \ij < fi'j; while for the region /zj > //J, the normal phase is unfavored and the pion 
superfluidity phase is the stable ground state. The distinctly different phase structure between 



14 



the chiral limit and the finite current quark mass suggests that the value of the critical isospin 
chemical has close relation with the pion mass. Even though for the simplicity of numerical 
study, the obtained critical point jjfj in this paper is not exactly the pion mass for the case 
beyond chiral limit, we point out that the improved calculation (more involved) within the 
DSE formalism should confirm jjFj = m w . 

Furthermore, our calculation shows that in the chiral limit with finite isospin chemical po- 
tential, both quark condensate and traditional mixed quark-gluon condensate vanish in OPE 
with the appearance of isospin density induced pion condensate and pion mixed quark-gluon 
condensate. In the real world, for m < fi c j, quark condensate and mixed quark-gluon con- 
densate exist in OPE with the vanishing of pion and pion mixed quark-gluon condensate; for 
fx j > /ij, the magnitudes of both isospin density induced condensates increase with the increas- 
ing of isospin chemical potential, while the magnitude of quark condensate decreases with the 
increasing of isospin chemical potential (It is expected that the mixed quark-gluon condensate 
has the similar behavior). Meanwhile, numerical calculations suggest that the induced pion 
condensate and pion mixed quark-gluon condensate have the same signs with their correspond- 
ing traditional chiral condensates. We also obtained that the magnitude of the ratio of pion 
mixed quark-gluon condensate to pion condensate is close to the one of traditional quark-gluon 
condensate to quark condensate in the vacuum, both for the chiral limit and for the real world 
in the pion superfluidity phase. 

Since there is no Fermion sign problem at finite isospin chemical potential with zero baryon 
chemical potential, in principle, the evaluation of the induced mixed quark-gluon condensate 
and four-quark condensate can be investigated through the lattice Monte Carlo method. The 
effect of these isospin chemical potential induced condensates on the hadron properties can be 
investigated in the framework of QCD sum rules. However, such effects on the hadron properties 
can more directly be explored using suitable Bethe-Salpeter equations in conjunction with the 
solutions of the quark gap equation. A natural extension of the present work is to investigate 
two flavor QCD phase diagram, condensates and hadron properties at finite isospin chemical 
potential for both nonzero temperature and baryon chemical potential. 
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